Where is the tt resonance? 
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We discuss the interplay of particle-particle and particle-hole spin-triplet channels in high- Tc su- 
perconductors using a quasiparticle dispersion motivated by angle- resolved photoemission. Within 
a generalized RPA, we find a well defined antibound state of two holes, the n resonance of Demler 
and Zhang, as well as a bound state of a particle and a hole, the spin exciton. We show that the 
energy of the tt resonance always exceeds 2A, twice the maximum d-wave gap, therefore the neutron 
resonance observed in the cuprates around energy A is most likely a spin exciton. At the same 
time, we speculate that the tt particle can exist at higher energies and might be observed in neutron 
scattering around 100 meV. 



74.20.-z, 74.72.-h, 76.50.+g 

Reflecting our poor understanding of higli-Tc super- 
conductivity in general, theoretical debates continue on 
virtually every aspect of it. A good example is the reso- 
nance observed in inelastic neutron scattering at energifg 
25-43 meV in YBazCuaOy-y and BizSraJpaCuaOg+aBEl 
Measurements using polarized neutronsH indicate that 
the observed excitation is magnetic in nature. Theo- 
retical proposals ifleJude a van Hove singularity in the 
Stoner continuum ,al3|a, hound state (spin exciton) in the 
particle-hole channel ,EI"EEI and an antitpimd state in the 
particle-particle channel (tt resonance) .lldO Deciding be- 
tween these approaches is somewhat difficult since in all 
theories the resonance mode is an excited pair of quasi- 
particles with totalJattice momentum Q = (tt/o, tt/a) 
and spin S = iMyta 

In this paper, we use available experimental data 
and a simple kinematic argument to demonstrate that 
a particle-hole bound state is the most likely explana- 
tion for the resonance. In a nutshell, our reasoning goes 
as follows. Angle-resolved photoemission spectroscopy 
(ARPES) directly measures the dispersion of fermionic 
quasiparticles ek in the superconducting state. From 
this, one can deduce information about the continuum 
of states of two quasiparticles with total momentum Q 
and, in particular, determine its energy bounds. An anal- 
ysis of ARPES data in Bi2Sr2CaCu208-i-5 shows that the 
two-quasiparticle continuum starts at E^i^ ~ 2A, where 
A is the maximum value of the superconducting gap at 
the Fermi surface. The commensurate neutron resonance 
in Bi2Sr2CaCu208+5 resides at £^res ~ A, definitely be- 
low the lower edge of the continuum. The resonance then 
can only be a bound state of two quasiparticles, not an 
antibound state. 

This argument — that the lower edge of the two-particle 
continuum is determined by A and not by the chemical 
potential /i as in Ref. |l^ — is based on the analysis of the 
Fermi surface inferred from ARPES data (Fig. ^. We 
have verified that, for this Fermi surface, E^i^ is deter- 
mined not by electrons along zone diagonals k^ = ky, but 
rather by electrons near the hot spots (points at k = \ip 



separated by Q). The latter are located close to (0,7r) 
and symmetry related points, and in the superconducting 
state have a gap A(khs) ~ A. This yields a threshold at 
£■111111 ~ 2 A. Note that this argument also invalidates the 
interpretation of the resonance as a-yan Hove singularity 
in the two-quasiparticle continuum.H 

Does it mean that the n resonance does not exist? Not 
necessarily. In our model calculations, the spin exciton 
and the tt resonance are found to coexist at intermediate 
coupling strengths, similar to earlier findings.E3 Further- 
more, in a situation where the continuum of states with 
two holes is narrow and has a sharp upper edge, we find 
that the tt resonance is rather sharp. To which extent 
these conditions are satisfied in the cuprates is not clear 
because self-energy effects could wash out the upper edge 
of the continuum. We therefore can only speculate that 
the TT resonance might be discovered in neutron scatter- 
ing at energies above E'min (60 meV in optimally doped 
Bi2Sr2CaCu208+5), more likely around 100 meV. 

To proceed, we need a model treating all spin-triplet 
excitations (with charges and ±2) on an equal foot- 
ing. The simplest approach is an adaptation of Ander- 
son's treatment of, phase fluctuations and plasmons in 
a superconductoitJ to the spin-triplet diajiael. Precisely 
this method, Jarown as generalized RPA,li3't3 was used by 
Demler et alS3. To make the approach self-consistent, we. 
make a quasiparticle approximation with a dispersiorJiZI 
fitted to peak positions in the ARPES data in the su- 
perconducting state of Bi2Sr2CaCu208-i-5. As the mea- 
sured dispersion is flat near (0,7r), our calculation yields 
a strikingly narrow two-electron continuum with a width 
of about 10 meV. This narrowness is caused by the prox- 
imity of the hot spots to the van Hove points at (tt, 0). 

This quasiparticle model is a "best case" scenario for 
the antibound state. This model assumes that one has 
a superfluid Fermi liquid, which appears to be consistent 
with transport and ARPES data, at least at low tem- 
peratures. Within this approximation, the upper edge 
of the two-hole continuum is sharp, which is favorable 
for the formation of the n resonance. By making such 



an approximation, we essentially neglect the large inco- 
herent part of the electron spectral function known to 
exist from photoemission. On the other hand, this in- 
coherent part is absent bRlaw-|the energy scale A -I- Ej-cs 
set by the spin resonance, I13'E3 so the fermionic incoher- 
ence affects the two-particle propagators only at energies 
above 2 A -I- -Ercs- Thus the tt resonance can, in princi- 
ple, live below this threshold. We make no attempt to 
explain why the quasiparticle dispersiao. is_so remarkably 
flat near (tt, 0). A number of authoralj CJ have pointed 
out that this flatness most likely results from a strongly 
w-dependent self-energy which renormalizes the bare dis- 
persion. Also note that the extreme narrowness of the 
two-particle continuum is not a crucial part of our argu- 
ment: we have observed similar behavior with a much 
wider (70 meV) two-particle continuum based on normal 
state dispersions. 

We now proceed with the calculations. Consider a sys- 
tem of fermions with dispersion, ej, given in Rcf. |lj and 
the nearest-neighbor interactionlid'Ell 



iv) 



VniTij. 



(1) 



The J term gives rise to an attraction in the particle- 
particle d-wave singlet channel and in the particle-hole 
triplet channel, where it can induce a spin exciton. The 
V term accounts for repulsion in the triplet pairing chan- 
nel and gives rise to the n resonance. In the conventional 
t-J model, V = 0. Following Refs. ^ and |l|, we con- 
sider a more general interaction and treat V and J as 
independent parameters. 

Although we are interested in spin susceptibility, in a 
superconducting ground state an operator of spin has the 
same quantum numbers as an operator creating a spin- 
triplet pair. We therefore|^alyze the linear response for 
the set of three operatorst3 
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7f = iV-V2^^^4 



Q-kT^-kT" 



(2) 
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Here a^u is an electron annihilation operator and g^ = 
cos{kxa) — cos{kya). Operators S+, n and n destroy 
a bosonic excitation with the same momentum Q = 
(7r/a,7r/a) and spin S^ — —1, but with different charges 
and ±2, respectively. 

As a warmup exercise, consider first a hypothetical 
(Fermi liquid) non-superconducting ground state. The 
triplet channels (0- 3) are decoupled by virtue of a charge 
U(l) symmetry. The Fourier transform of the bare pair 
susceptibility 



X^(t) = -z0(t)([7r(t),^t(o)])| 



y=j=o 



(5) 



is shown in Fig. |2|(a). As anticipated, the two- hole con- 
tinuum [uj < 0) is strikingly narrow: its upper edge is 
at |i?max| "^ 10 meV. The uj^^/'^ divergence near E^ax 
is a van Hove-type singularity associated with the (tt, 0) 
points. Observe that, at |a;| > |i?inax|, X^i^) is purely 
real. On the other hand, the bare spin susceptibility 
X%i^) [shown in the insert of Fig. 0(a)] is complex for all 
frequencies. 

Within the generalized RPA approximation, the full 
susceptibilities are given by 



XTr(w) - 
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X°(^) 
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(6) 



Note the opposite signs of the interaction terms: attrac- 
tive in the particle-hole channel, repulsive in the particle- 
particle channels. Because x'ii^) = V^ ^^^ a solution 
for any V > [Fig. ||(a)], the full n susceptibility X7r(i^) 
acquires a pole above the upper edge of the continuum. 
This pole is the tt resonance. In contrast, no pole oc- 
curs in the RPA response for 5"+. This is a consequence 
of the fact that the particle-hole continuum [and, hence, 
Imx5(a;)] extends to the lowest energies. 

We now turn to the superconducting state. To fa- 
cilitate the RPA treatment, we assume that the super- 
conducting state is of the BCS type with a d-wave gap 
A(k) = A{cos kxa — coskya)/2. We treat A as another 
input parameter in the problem. In principle, the gap 
has to be computed self-consistently, but for our pur- 
poses this is not necessary as we are not concerned with 
response functions in the singlet channel at small mo- 
menta. 

Since in a superconductor charge is defined modulo 
2, all three operators (0-0) now carry identical quantum 
numbers and one can use any superposition of these. It 
is customary to choose 



5+, Ai 



A^ 



V2 ' ^" V2 
The operators Ai and A2 describe fluctuations of t 



(7) 



phase and the amplitude of the tt mode, respectively.E3 
The bare susceptibilities xj? , {u>) are defined, similarly to 
Eq. (P), as the Fourier transforms of 

x',,,(t) = -t0{t){[A,,{t),Al{O)]) 



iv=j=o 



(8) 



They can be written as components of a 3 x 3 matrix x'^. 
At T = 0, we have 



C'H = ^e( 




where E^ = e^ + Aj., k' = Q — k, and 9!), are components 
of the vector 

0= (2-i/2p(k,k'), 5kri(k,k'), 5k^(k,k')) , (10) 



where I, n and p are the BCS coherence factorsE 

Z(k,k') = MkMk' +WkWk', 
n(k, k') = Wk^k' - WkWk', 
p(k, k') = UkWk' + WkWk'- 



(11) 



The momentum sums are evalulated directly by including 
a small broadening w — > w + zF in the energy denomina- 
tors. Note that, due to strong particle- hole asymmetry 
implied by our dispersion, ofF-diagonal terms in x^ are 
by no means small. 

Within the generalized RPA scheme, the full suscepti- 
bilities Xqq' (^) ^re given by the matrix RPA equation 



xH = x"M + xH^rH 



(12) 



where V = diag(— 2J, V, V) in the basis (^. Recall that 
both J and V are assumed to be positive. We next dis- 
cuss the behavior of response functions in various situa- 
tions. 

J j^ 0, V rffirfh This is the situation considered in 
most studies. qBO'lI Although the BCS condensate now 
mixes all components of Xqo'(^)i the full spin suscepti- 
bility xoo('^) = X£,('Vi) ^^ ^till given by the standard RPA 
expression, Eq. (^).lI This effective decoupling is due to 
the fact that the J term in ™) yields no interaction in the 
particle-particle triplet channel. Despite a formal simi- 
larity with the normal state result, the form of Xoo('-^) 
changes dramatically in a d-wave superconducting state: 
the Stoner continuum develops a hard gap 

-Bn,in(Q) = min(£k + £^Q-k) = 2|A(kh,)| « 2A. (13) 

The imaginary part of Xool'^) jumps from to a finite 
value at 2A(khs). By Kramers-Kronig relation, this dis- 
continuity causes a logarithmical divergence of Rexoo(^)- 
Therefore, for adpittary J, the full xoo('-^) h^^s a resonance 
pole at w < 2AB113 

In the absence of the triplet-pair coupling V, there is no 
TT resonance above the continuum. At the same time, the 
TT-TT correlation function contains a resonance below the 
continuum, which is nothing lmt_the spin exciton mixing 
into all other triplet channels.&Ej Indeed, for V — 0, the 
solution of Eq. (n3) can be written in the following form: 

xqA^) = x°A^) - 2x^o(^) J'^^'H^) x!^,'(^), (14) 

where J^^^{lu) = J/[l + 2Jxoo(^)]- The spin-exciton 
pole in the RPA vertex J^^^{uj) shows up in all full sus- 
ceptibilities. 

J = 0, V y^ 0. In this limit, there is no bound state in 
the particle-hole {q = 0) channel because J = 0. The tt 
resonance in the particle-particle channel, which already 
existed in the normal state, is pushed to a higher energy 
when the superconducting gap opens up. We stress once 
again that the lower boundary of the two-hole continuum 
^min is produced by fermions with momenta near hot 
spots, hence i?min ~ 2A. Certainly, an antibound state 



in the particle-particle triplet channel is located above 
2A. How much above 2A depends on the width of the 
two-electron continuum and the coupling strength. 

The bare tt-tt response function in a superconducting 
state is presented in Fig. 0(b). It has two step- like dis- 
continuities at -Emin and -Emax, which are seen in both 
hole-hole (w < 0) and particle-particle {ui > 0) spectra. 
The location of the tt resonance in a superconductor is, 
however, not simply given by (ph, but is a solution of the 
secular equation 



det 



xU^) x^2(^) 



= 0. 



(15) 



|i 



All four matrix elements x^q' for q = 1,2 have the same 
thresholds as xItt i^ Fig g(b) . We solved this equation 
numerically and found that the tt resonance indeed moves 
to an energy higher than 2A. For V — 40 meV, the 
resonance moves from 22 meV in the normal state [Fig. 
^(a)] to 82 meV in the superconducting state with 2A = 
70 meV [Fig. g(b)]. Reasoning similar to that above 
implies that the tt resonance shows up in the spin channel. 

J y^ 0, V ^ 0. The general case interpolates between 
the two limits. Fig. H presents our results for the spin 
susceptibility at various ratios of the coupling strengths 
V/J. Interestingly enough, for comparable couplings, we 
found that the response functions contain two peaks at 
different energies. One, above E-a^ay,, is the tt resonance, 
the other, below -Emin, is the bound state in the particle- 
hole channel. As V is increased, the high-energy reso- 
nance pulls away from the upper edge of the two-hole 
continuum and stregthens. At the same time, the low- 
energy resonance approaches the lower edge and enters 
into the two-electron continuum, where it broadens and 
finally becomes invisible. In the opposite limit, as V gets 
progressively smaller, the tt resonance weakens, merges 
with the outer edge of the continuum-and disappears. 

Discussion. Direct comparisonE3 of neutron and 
ARPES data for Bi2Sr2CaCu208+5 shows that the neu- 
tron resonance is located well below 2A. We believe 
that the same holds true for YBa2Cu307_y. Because 
the lower continuum edge £^min is just below 2A, the 
resonance almost certainly occurs below this edge. This 
point can be verified by analysing ARPES data at the hot 
spots. In a slightly overdoped Bi2Sr2CaCu208+i5 sample 
with X. = 87 K, the gap value at the hot spots is 32 ± 3 
meV,cJ and Ey^m > 58 meV. The neutron resonance at 
Q = (tt/c, 7r/a) in the overdoped Bi2Sr2CaCu20*|jT5-|With 
similar T^ = 83 K is observed at E = 38 meV,MEl i.e., 
well below E^min- We thus conclude that the neutron res- 
onance is an exciton-likc bound state in the particle-hole 
channel, not an antibound state, such as the tt resonance. 

It has been previously remarked that particle-hole.piipL=. 
ing can mask the origin of the neutron resonance all3'Lj 
a spin resonance mixes into triplet-pair channels and 
vice versa. We rely on a continuity argument to make 
the unambiguous identification of the experimental neu- 
tron peak, which is observed below the continuum edge 



£"111111 ~ 2 A, as expected for V <^ J. As the pair-triplet 
coupling V is turned off, the low-energy peak stays below 
the edge and continuously evolves into the usual RPA 
spin resonance at V = 0. On the other hand, there is 
no collective mode below £^min in the limit V ^ J, in 
which the n resonance can be defined unambiguously. 
Moving the high-energy tt resonance across the contin- 
uum involves a discontinuous change. For this reason, 
the neutron resonance should not be associated with the 
TT resonance. 

We feel that a more realistic analysis will not change 
this general conclusion because our argument does not 
rely on detailed dynamics of the resonances, but is based 
on measured kinematics of low-energy fermion excita- 
tions and on general properties of bound states. If one 
wishes to find the tt resonance in neutron scattering, the 
search should be confined to energies above the two-hole 
continuum. As its true upper edge is not known exper- 
imentally, we can only suggest that this energy exceeds 
2 A. On the other hand, because of the large incoherent 
parts of the spectral function observed in ARPES data 
for energies beyond. A, there may be no true upper edge 
to the continuum,E2l and it is quite possible that the it res- 
onance will be strongly damped, or perhaps even absent, 
in neutron data. 
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FIG. 1. "Hot spots" are points on the Fermi surface con- 
nected by the wave vector Q — (vr/a, n/a) or those equiv- 
alent to it. Equivalently, they are given by intersections of 
the Fermi surface with the magnetic zone boundary (dashed 
lines) . 
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FIG. 3. Spin susceptibility [— Inixs(a')] in the supercon- 
ducting state at Q. Values of V are shown near each curve. 
-Bmin and -Emax are boundaries of the two-hole continuum. 
The chosen values for J and A are similar to those measured 
by neutron and ARPES, respectively. 
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FIG. 2. The bare response function (T = 0) of the operator 
TT in a hypothetical Fermi liquid (a) and a d-wave supercon- 
ductor with Amax=35 meV (b) at Q. Solid line: real part, 
dashed line: imaginary part. Insert: spin susceptibility in the 
normal state (note the difference in vertical scales). A broad- 
ening of r=0.5 meV was employed. Solution of the equa- 
tion x^i'-^) ~ 1/^ gives the RPA energy of the tt resonance, 
oj = —Ett. Note the difference in vertical scales. Horizontal 
arrows indicate two-particle continua. 
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